Possible spin triplet superconductivity in Na^CoC^ ■ yH 2 
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Combining symmetry based considerations with inputs from available experimental results, we 
make the case that a novel spin-triplet superconductivity triggered by antiferromagnetic fluctuations 
may be realized in the newly discovered layered cobaltide Na^OoC^ ■ J/H2O. In the proposed picture, 
unaccessable via resonating-valence-bond physics extrapolated from half-filling, the pairing process 
is similar to that advanced for S^RuCU, but enjoys a further advantage coming from the hexagonal 
structure of the Fermi-surface which gives a stronger pairing tendency. 



Anderson's suggestion that a doped antiferromagnet 
on a 2d square lattice be best described in terms of res- 
onating valence bond (RVB) states [1], set the direc- 
tion for a huge subsequent effort pursuing this picture 
in an attempt to understand superconductivity in the 
cuprate oxide compounds. While this has lead to the 
development of new avenues in the physics of low di- 
mensional condensed matters, conclusive evidence of such 
novel states has not turned up. Meanwhile, renewed in- 
terests in RVB physics in systems on triangular lattices 
(TL) have been emerging recently - though the S=l/2 
antiferromagnetic (AF) Heisenberg model on a TL is gen- 
erally regarded as having a ground state with long range 
order, closely related systems such as the undoped quan- 
tum dimer model [2] and the t — J model on a TL with 
the hopping sign t > [3] have been shown to exhibit a 
certain amount of stability of RVB-like states. 

The discovery by Takada et al of superconductivity in 
Naa;Co02 -y H2O (x ~ 0.35, y ~ 1.3) [4] , a system made 
of stacked layers of cobalt ions sitting on a 2d TL and 
surrounded by octahedrally displaced oxygen atoms will 
no doubt spark new interest in this line of investigation. 
Superconductivity was achieved at T C =5K by oxidation 
of the sodium ions (changing the effective charge den- 
sity) followed by intercalation of the C0O2 layers in the 
mother compound NaCo 2 04 with H 2 0, thus effectively 
enhancing the two dimensionality of the system. Guided 
by analogy with the cuprates, the authors of ref [4] were 
lead to suggest that this system maybe understood in 
terms of doped AFs on TLs, which may well be related 
to RVB states. 

As much as such directions are tempting, one must 
proceed with due caution since the experimental situa- 
tion is rather complex. As discussed shortly, the sys- 
tem is actually away from the vicinity of half-filling, the 
regime where RVB physics is considered to be valid. It 
is reasonable then to resort to generic methods no more 
specialized than a combination of symmetry analysis and 
simple "fermiology" , reinforced by experimental inputs. 
While experiments conducted on Na x Co 2 -y H 2 are at 
present basically limited to resistivity and preliminary 
susceptibility data, we have noticed that several conclu- 
sions -which are in fact remote from what an RVB-based 



investigation would lead to- can already be drawn. We 
therefore believe it useful at this stage to raise notice on 
this point with the hope of motivating further activity. 
The seminal work of Rice and Sigrist, who adopted a 
similar approach in the case of Sr 2 RuC>4 [5], had played 
an instrumental role in the subsequent identification of 
the triplet superconductivity in that compound. 

During the final stage of this work, we came across a 
preprint by Baskaran [6] addressing the same system. Al- 
though his reduction of the problem to a one-band model 
is similar, the discussion there is based on the RVB pic- 
ture. The present work is founded on a complimentary 
approach as explained above. 
Deduction of effective theory 

The first point to make is that the present system can- 
not be taken as a realization of the dilutely doped AF 
on a TL. The 3d cobalt ions in the parent compound 
NaCo 2 C>4 with a formal valence of +3.5, are actually a 
mixture of the low-spin states Co 3+ (t 2g eg with S=0) and 
Co 4+ (t 2g e g with S=l/2). With no experimental indica- 
tions of charge ordering in NaCo 2 04, the most natural 
starting point for studying these systems is as a random 
mixture (1:1) of S=0 and S=l/2 states distributed on 
the TL. Hence instead of starting near half- filling, we 
must consider the case of three-quarter filling, which is a 
totally new situation. Superconductivity occurs after ox- 
idation of the sodium ion so that the compound becomes 
Nao.35Co0 2 -y H 2 0, in which case the system is approx- 
imately 1/3 filled. We therfore stress that RVB physics, 
presumably valid close to half-filling, by no means has a 
trivial extrapolation to this regime. 

The next issue to consider is the magnetic properties 
of the system. The susceptibility [7] of NaCo 2 C>4 con- 
tains a Curie- Weiss component with large negative Curie 
temperature Tc U rie ~ —285K, a trend also confirmed 
in measurements on Nao.3sCo0 2 -y H 2 0, though with a 
reduced value of Tcurie ~ —18.8K [8]. Taken together 
with the finding of spin density wave (SDW) like ten- 
dencies upon Co to Cu subsitution in NaCo 2 04 [9], it 
is natural to assume a fair amount of AF fluctuation to 
be present. LDA calculations show that the Fermi sur- 
face for NaCo 2 04 sits near the top of the ai s like state, 
split off from the t 2s manifold due to oxygen distortions 
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[10]. These circumstances suggest a one-band Hubbard 
model (U w 5 — 8eV, t leV) on a TL around quarter to 
third filling as the simplest model capturing the essential 
physics of this system. At this point one may already 
expect nontrivial pairing tendencies. For instance, there 
are perturbative studies that claim that triplet pairing is 
favored over singlet pairing [11]. With direct quantum 
Monte-Carlo methods unavailable due to severe negative 
signs, we proceed with the generic approach mentioned 
earlier. 

Symmetry considerations 

Group theoretical studies, valid irrespective of the pair- 
ing mechanism, have played essential roles in classifying 
superconductivity in heavy fermions, organics, and e.g. 
the ruthenate compound [12,13]. For the present case, 
the relevant symmetry group is Q = Dq^ <g> U(l) ® T, 
with D^h the hexagonal group together with the inver- 
sion symmetry about the plane, U{1) the gauge symme- 
try broken by the onset of superconductivity, and T the 
time reversal symmetry. For Nao.3sCo02 -y H 2 0, mag- 
netization measurements suggest the presence of a mag- 
netic anisotropy [8]. Hence the absence of SU(2) from 
Q. The irreducible representation for this situation can 
be classified into 12 classes, r„,l < a < 6, with +(- 
) standing for spin singlet (triplet) states. States with 
subscripts 1 < a < 4, and those with a = 5, 6 are Id and 
2d representations respectively. The conventional point 
of view [14] that AF fluctuations are most compatible 
with singlet pairing will lead us to one of these (+) rep- 
resentations. 

A more intriguing possibility however is the case of 
triplet pairing. Not easily accessed from RVB treatments 
which deal primarily with spin singlets and their fluctu- 
ations, it is perhaps worthwhile to highlight their exotic 
properties. Later on we will mention that experiments so 
far seem to spell out such a state as the one most prob- 
able, and a novel scenario from which we recover this 
pairing choice will be described. As mentioned there, we 
will be interested in states with the d-vector perpendicu- 
lar to the cobalt plane, i.e. dj / z, which ensures that the 
electron spins contributing to the wavefunctions have a 
vanishing z-component. It turns out that the repre- 
sentation realizes such states. The fourth order terms of 
a Ginzburg Landau expansion of the free energy is of the 
form 

T q uad = J dM/MM 2 + \V2\ 2 f + - VlV2)% 

(1) 

where 7/1 and 772 are coefficients of the two basis func- 
tions zk x and zk y . When the coefficient j3 2 > 0, we have 
d = z(k x + iky), a state violating T. The same pair- 
ing state also arises in the tetragonal (D^h) case, though 
additional conditions on the GL coefficients are required 
there. The T-violation is most dramatically manifested 



in electric/magnetic properties. The derivative terms of 
the GL theory for representations Tf are 

T = J dr\K x \D x n Y + D yV2 \ 2 + K 2 \D x v 2 - D ym \ 2 

+K 3 (\D xVl - D yV2 \ 2 + \D xV2 + D ym \ 2 )} (2) 

with = — 2ieA ll the covariant derivatives. 

All z-derivatives are omitted on account of the two- 
dimensionality. Adapting a method due to Furusaki et al 
[15] who discussed the tetragonal case D^, we find that 
the terms in the GL Lagrangian density coupling linearly 
to the electric field E read 

C E = c\{E x r){ + E y ii 2 )(D x rii + D y n 2 ) + c.c. 
+ c 2 (E x r] 2 - E y nl)(D x n 2 + DyTfi) + c.c. 
+ c 3 {(E x nl - E y r)2){D x r]i - D y r] 2 ) 

+ {E x rf 2 + E y nl){D x n 2 + D ym ) + c.c.} (3) 

where the coefficients Cj (l<j<3) are proportional to the 
Kj's in the free energy. From these one extracts the fol- 
lowing contribution; 

£cs-iikc = *( c 3 - ci - c 2 ){nln 2 - rj^Vi) 

x(A y d x A - A x d y A ), (4) 

which is of the Chern-Simons (CS) form. Several con- 
sequences follow, which are basically similar to what 
has been discussed in the context of SrRu 2 04. First, 
a zero field Hall effect is expected, due to a sponta- 
neous (orbital) magnetization d ■ (k x V^)rf cx if* x ff , 
where ff = (771,772)- Futhermore one sees from the 
nonunivcrsality of the coefficient of the CS term that 
the Hall conductivity a xy is not quantized. The phys- 
ical reason for the latter feature can be identified micro- 
scopically. (Readers are referred to refs. [16,17] for the 
tetragonal case.) Inserting the d-vector d = Ao(/i(fc) + 
if 2 (k))z, with fi(k) = 2 sin(^k x ) cos(|fc y ) and f 2 (k) = 

2cos(—^Y'k x )sm(^ky) + 2sinfc y into the Bogoluibov-dc 
Gcnnes Hamiltonian 

H = ^[(e(k) - M )t 3 ® 1 + n ® i(d{k) ■ f)r 2 ]*, (5) 

where e(k) = —4 cos(^-k x ) cos(|fc y ) — 2 cos(fc y ) and per- 
forming a gradient expansion to one-loop order, one sees 
that the correction to the universal contribution o xy = 
e /4-7T comes from the normalization of the generating 
functional due to the Goldstone mode associated with the 
breaking of the gauge symmetry. This simply reflects the 
fact that the charge is not a well-defined quantum num- 
ber in a BCS condensate [17]. On the other hand, the 
residual U(l) degree of freedom left of the SU(2) spin ro- 
tational symmetry suggests that the transport of S z has 
universal feautures. This expectation is verified by gaug- 
ing the system via the coupling to a spin-gauge field, from 
which one confirms that the spin Hall conductivity a s 
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is quantised in integer multiples of l/27r7i. The integer 
factor is a topological invariant related to the chirality of 
the superconducting state, characterized in terms of the 
quantity g = '(Rerf z ri, — Imrf z ri, el) as 



N, 



chiral 



d 2 k trjgj V fcx g x V fca g) 
(2tt) 2 t r|g|i 



(6) 



which is +1 for the present choice of d-vector (-1 for the 
opposite chirality). To complete the discussion on sym- 
metries and exotic topological properties, we note that 
vortices obeying nonabelian statistics can be realized in 
this chiral p-wave state [18]. 
Relation to experimental inputs 

We list the salient features of the inputs provided by 
experiments conducted so far [8]. 

1. The normal state resistivity as a function of tem- 
perature fits rather well to a T 2 -dependence. 

2. Preliminary NMR measurement of T^ 1 at the 
cobalt site shows a prominent peak just below T c , 
reminiscent of the Hebcl-Schlichter peak for con- 
ventional s-wave superconductors. 

3. Spins direct their moment in-plane under an ap- 
plied magnetic field ~ 8T. 

Deferring discussions on feature 3 for the moment, we 
focus on the first two. We take feature 1 to be an indi- 
cation that the system lies in the weak coupling regime 
[5]. It then becomes plausible to interpret feature 2 in 
terms of a fully-gapped pairing state without nodes. Al- 
though this is usually considered a hallmark of s-wave 
pairing, it is worth recalling that the Balian-Werthamcr 
state, with oIbw — A fc has a sufficiently strong singular- 
ity of the density of states to exhibit a coherence peak as 
well. The chiral p-wave state (with the cylindrical Fermi 
surface) discussed above is the two-dimensional analogue 
of this situation. An in-plane measurement of the phase- 
sensitive thermal conductivity on single crystals, when 
they become available, is expected to give an isotropic 
result. (The absence of a peak in the case of Sr2Ru04, 
which had also been considered to be a chiral-p wave 
superconductor, is consistent with the existence of hor- 
izontal nodes related to interlayer pairing [19].) States 
such as d — p x z and d = p y z, permissable on symme- 
try grounds, can be dismissed on bases of these observa- 
tions. Alternative fully-gapped states include the singlet 
d x 2_ y 2 +id xy and d x i_ y i +is states [6]. One can directly 
distinguish between singlet and triplet states by perform- 
ing a Knight-shift measurement. For singlet pairings, the 
shift K will drop to zero on approaching T = for all 
magnetic field directions, with K cx Y(k,T) w e~~ , 
where Y(k, T) is the Yoshida function. Triplets on the 
other hand will exhibit an anisotropic magnetic suscep- 
tibility in the superconducting state; an in-plane field 



ffli will polarize both the Cooper pairs and the quasi- 
particles, so that K//(T) = const, is expected, whereas 
the shift for H//z will lead to a behavior similar to 
the singlet case, K± cx Y(k, T). Remarkably, prelimi- 
nary Knight shift measurements do show a temperature- 
independent behavior [8]. (A large H c2 ~ 61T [8] also is 
suggestive of triplet pairing.) These results, taken at face 
value single out the nodeless triplet pairing, of which the 
chiral-p state is the simplest. 

Possible scenario for triplet pairing in Naa.^Co02 -y 
H 2 

According to LDA calculations [10], the central cylin- 
drical Fermi surface of NaCo 2 04, with dominant a\ g 
character bears a shape close to a hexagon when viewed 
in the k z = const, plane. Hence there are three (approx- 
imate) nesting vectors, given by Qi = ^TKi,i = 1,2,3. 
Here we define the if-points as Ki — C^ 1 {\G\ + 3G2), 



Gi = 2tt 



bxc 



G 2 
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in which a = (1, 0, 0), b 



a-bxc ' ^ 2 a-bxc ' 

(— 1/2, -s/3/2, 0), c = (0,0,1) are the primitive vectors of 
the layered TL. Cq denotes the anticlockwise rotation 
about the z-axis by an angle of n/6. The presence of 
such nesting instabilities in the T — K direction has been 
suggested by Terasaki et al [9] based on their observa- 
tion of an order-from-disorder type emergence of a spin 
density wave (SDW) in Cu-substituted NaC01.gCuo.2O4. 
We are assuming below that these qualitative Fermi sur- 
face features remain intact for Nao.3sCo02 -y H 2 0, which 
appears to be in accord with the observed Curie- Weiss 
contribution to the susceptibility. 

Fermiology, the study of many-body effects with an 
emphasis on Fermi surface properties has been playing 
a central role in advancing our knowledge on supercon- 
ductivity in various correlated electron systems [20] . One 
such scenario [21], which takes advantage of nesting ten- 
dencies and leading to triplet pairing, has been advanced 
for the tetragonal lattice case, with the aim of explain- 
ing the pairing in Sr 2 Ru04. We would like to point 
out that a very natural extension to the present TL 
case is possible. This argument starts by recasting of 
the linearized BCS gap equation at T c into the form, 



T c (x e jv <°)<< v «>>fs j w here we are following the nota- 
tions of ref. [22]: << V$ »fs is the pairing interaction 
V(k — k') averaged over the Fermi surface, i.e. 



« Va, »fs- 



J FS dk J FS dk'V^k- k')(t>(k)<i>(k') 



[J FS dk']J FS dk^(k) 



(7) 



with 4>(k) the /^-dependent part of the order parame- 
ter, i.e. A(fc) = A <p(k). We may view this as a 
sort of variational principle in which the actual pairing 
occurs for cases which yield a (l)positive and (2)large 
<< V,j, »fs- In the present situation, where we 
are dealing with AF fluctuations, the pairing potential 
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V(Qi) > 0. For this to lead to pairing requires that 
4>{k)4>{k + Qj) < for some i. The simplest way to real- 
ize this is to introduce a nodal direction which coincides 
e.g. with the k x axis, i.e. a p y -pairing state. One im- 
mediately sees that this choice is highly beneficial as all 
three nesting vectors satisfy the required condition (see 
Fig. 1(a)), and therefore the entire Fermi surface is avail- 
able for the pairing. (For this Pj,-pairing, the Qi-channel 
should be the most dominant because it coincides with 
the direction in which the magnitude of the gap becomes 
maximal.) In this respect, the advantage of such node 
formation (vertical nodes on the central Fermi surface) 
is bigger than in the ruthcnatcs where only one pair of 
nested segments is involved in the pairing. This state is 
degenerate (i.e. has the same T c ) with those in which 
the nodes are located at the intersections between the 
hexagon Fermi surface and the lines k y — ±\/3k x (Fig. 
1(b) and (c)). Further gain in condensation energy is 
achieved by constructing a linear combination of these 
three states in such a way that the nodes will cancel one 
other, opening a full gap at all points on the Fermi sur- 
face. This becomes possible through the construction 



0)p y +ip* cx <\> Vv (k) + i (j) Py (C 6 k) + <\> Vv (C 6 2 £)) . (8) 




FIG. 1. Illustration of pairing via nesting of Fermi surface 
segments, (a) Node located at p y —0, corresponding to p y pair- 
ing, (b) and (c) denote the case where the nodal directions 
are rotated by 60 and 120 degrees respectively. Superposition 
of the three states, with relative phases of \ between (a) and 
(b), (c) will cancel the nodes, leading to a chiral p-wave state 
with enhanced condensation energy. 

As will be discussed elsewhere, one can check that a 
similar construction cannot be made for an f-wave state. 
The appearance of the chiral p-wave state above relied 
on the Fermi surface geometry (reflecting the hexagonal 
symmetry), and the assumption of a nesting tendency. 
It is desirable that fine quality samples, amenable to de- 
tailed neutron scattering measurements be made in or- 
der to make a direct verification on this point possible. 
Further support for triplet pairing comes from the afor- 
mentioned magnetic anisotropy; when we interpert this 
in terms of an easy-axis along the c-axis, evaluation along 
the lines of ref. [22] can be used to demonstrate that the 



anisotropy ratio of the in-plane and longitudinal suscep- 
tibilities will directly enter into expressions for the triplet 
and singlet pairing potentials, and a sufficient anisotropy 
will eventually favor triplets (with dj ' / z) over singlets. 
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